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ABSTRACT 


The sensitivity of the coupled chemistry/diffusion model's 
outputs to a wide range of variation of the model's Independent 
variables has been investigated. It is shown that the efficiency 
with which the NO^ catalytic cycle destroys ambient Oo is extremely 
sensitive to the amount of NO emitted and to the relative rates of 
turbulent diffusion and chemical reactions. For representative 
conditions in the stratosphere, a tenfold variation of either the 
turbulence intensity or the reaction rate constant or the source 
strength can vary the efficiency from \% to 50^, If the duration 
of Phase III is a significant fraction of the total residence time 
of the plume, then these efficiency variations can alter Oo deple- 
tion rates by more than a factor of two. These results, therefore, 
point toward those variables which must be accurately defined or 
measured if one is to adequately predict the effect of SST opera- 
tions on the ambient Inventory of 0^ in the lower stratosphere. 
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INTRODUCTION 


Since the summer of 1972, A.R.A.P. has developed a second- 
order closure scheme for Isothermal reactions in a binary chemical 
system. This chemical kinetics model has been successfully coupled 
with the Invariant turbulence and diffusion models for the case of 
two-dimensional diffusion of two reacting chemical species which 
emanate from common or separate sources. The coupled model, while 
still necessarily an approximation, retains the second-order corre- 
lations for both species concentrations and turbulent motions and 
therefore permits separate evaluation of the flux divergences and 
inhomogeneous chemical reaction rates In the mass conservation 
equations for the reacting species. References 1 and 2 describe 
the development of this coupled model and the, results of some 
preliminary applications to the SST exhaust plume. We will try to 
avoid repetitions of those discussions In this report. 

The coupled model represents a substantial improvement In the 
simulation of chemical reactions in turbulent flow fields. Conven- 
tional techniques may fall to adequately describe the fate of SST 
exhaust products because of two severe limitations. First, because 
these techniques employ first-order closure and assume that the 
local reaction rate proceeds according to the local mean value of 
the species concentrations, they fall to account for spatial Inhomo- 
geneitles of the two reactants and the possible reduction of the 


r eacti on rate due to the correlated concentration fluctuations 
(C^Ci). In contrast, second-order closure permits dynamic simula- 
tion of this "mixedness" term and thereby permits quantitative 
evaluation of the effects of spatial inhomogeneity. Some of the 
work reported here demonstrates a mixedness limitation approaching 
25^ for equilibrium conditions and even larger effects for non- 
equilibrium conditions and in the core of the plume. Whether such 
an effect is significant depends on the uncertainty one wishes to 
attach to each of the indepe'ndent -variables 'of the model. 

A second and more serious limitation of conventional tech- 
niques for stratospheric plumes is the decoupling of the chemical 
depletion and turbulent diffusion processes and the assumption 
that the reactants are always intimately mixed. This assumption 
completely ignores the experimental evidence that the stratosphere 
is stably stratified and, for the most part, is in nearly laminar 
motion. It is quite possible that due to very slow turbulent 
diffusion the ambient species may be rapidly depleted by chemical 
reaction in the center of the plume and only slowly replenished by 
diffusion. As a result, significant portions of the plume may be 
shielded from the environmental species (as long as plume concen- 
trations remain high) and such plumes may survive considerably 
longer than previously estimated. 

The second-order closure model may be utilized in situations 
ranging from those in which vigorous mixing and slow chemistry 
make the flux divergence completely dominant, to situations of slow 
diffusion and fast chemical reactions in which the reaction rate 
dominates. Preliminary results, simulating the reaction of Oo with 
NO in the post -vortex SST exhaust plume have demonstrated a strong 
diffusion limitation on this reaction in the stratosphere - a 
condition which may persist for days to weeks. 

The work reported upon here is an extension of these early 
results in the form of a sensitivity analysis of the model's 
response to systematic variations of the input variables. In 
addition to the determination of the sensitivity of the model out- 
puts to the input parameters, and therefore knowledge of how well 
these input parameters must be known for reliable results, the 
choices of the ranges and values of each input parameter based on 
realistic values provides a first "catalog" of the expected ranges 
and interactions of diffusive and chemical processes in the SST 
far wake. These results may both guide and be verified by the 
NASA stratospheric jet wake experiments. They will also provide 
guidelines as to when conventional decoupled chemistry and diffu- 
sion models may be used and when coupled models are required for 
useful estimates of the far wake behavior. Finally, these results 
will be extremely useful in structuring three-dimensional models. 
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SYMBOLS 


E 

E 


f . 

,1 

k 

K 

P 

q 

Q 


r' 

* 


R 


R1 

S. 


d 


R 


u 

V 


constant evaluated using Eq. ( l8) 
concentration of 1th chemical species 
chemical efficiency defined by Eq. (33) 
average efficiency for residence time of plume 
transfer function 
reaction rate constant 
constant .In 'Eq, (6) 

diffusion limitation defined by Eq, (34) 
mixedness limitation defined by Eq. (35) 
any dependent variable 

twice the square root of the turbulent kinetic energy 

amount of NO emitted per unit distance of travel and 
per unit width of plume 

total reaction rate in plume defined by Eq. (31) 

maximum reaction rate in plume 

Richardson number 

critical Richardson number 

source or sink of a chemical species 

slope of iso-x^^^ contour in nomogram 

time • 

time required to dilute plume by 10-^ 

potential temperature 

residence time of plume 

longitudinal velocity 

lateral velocity 
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w 


X 

y 

z 


k 





V 

a 


vertical velocity 
longitudinal distance 
lateral distance 
vertical distance 
Arrhenius activation energy 
pDltzmahnl.cohstant- • 
diffusive scale length 
turbulence scale length 
kinematic viscosity of fluid 
standard deviation 


superscripts 

— mean component 

' fluctuating component 

subscripts 


chem 

e 

a 

P 

y, 6 
M 


depletion rate due to chemical reactions 
environmental value 
0^ chemical species 
NO chemical species 

reactive chemical species other than or NO 
unreactlve chemical species 


THE COUPLED DIFFUSION/CHEMISTRY MODEL 


The coupled model was developed for the following conditions 

1. Incompressible flow 

2. Binary and Isothermal chemical reactions 

3. Horizontally homogeneous turbulence 
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4 . 


No coupling between the chemical reactions and 
the turbulence field. 

5. Horizontal straight-line flow with the mean wind 
speed, U , a function of height only. 

6, Two-dimensional (z and t;) coupled diffusion and 
chemistry (l.e. lateral and longitudinal diffusion 
negligible in comparison with the vertical component). 

Despite these apparent restrictions, the model does capture a 
large part of the SST wake features and does permit the first 
critical study of the interaction between diffusion and chemistry 
in a turbulent, inhomogeneous system. In the future, each of these 
constraints will be removed, although the additional computational 
complexity will be considerable. It therefore behooves us to 
exercise the present model vigorously in order to pave the way for 
this future work. 

We consider the irreversible binary chemical reaction between 
chemical species a and p to form y and 6 . 


a + P — ► 7 + 6 


( 1 ) 


The constituents satisfy the species conservation equation 




(2) 


where C. is the concentration of the 1th chemical species and M 
represents all nonreactants in the reaction volume. We assume that 
the reaction rate for any joint value of the concentrations is given 
by * 


’ac 

a 

at 


chem 


= - VaS 


(3) 




chem 




and that 


(4) 


a 





(5) 


where kj^ are the reaction rate constants and are given by the 
Arrhenius equation 
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( 6 ) 


■ = KT^ exp (eJ/cT) 

Equations (3) and (4) represent the local instantaneous rate 
of change of the concentration due to chemical reaction. In order 
to determine the average rate of change, we assume the history of 
the joint values of and Cg at a fixed location comprises a 

stationary time series and each may be dissected into Its mean and 
fluctuating components 

0, = 0^ + C. ^ (7) 

and, by definition. Cl = 0 . Under these assumptions, the equa- 
tions for the average rates of change of a and p become 


'ac 


a 




chem 


= -k, (C Cft + CCA) 
1^ a p a p' 


( 8 ) 


ac 


p 


dt 




-I chem 

and similarly the species conservation Eq. (2) becomes 


C ' 
a 


(9) 


+ Cp + C^ C^ + C^ = 1 

(10) 

^ C + c^ 1 C-J = 0 

(11) 


Now, the equation for the conservation of mass of the a species 
may be written 

DC ac ac ac ac 

a a a a a 


a^c 


= V 


a 


a^c 


a 




a 


ax' 


8 y' 


az‘ 


ac 


a 


-“chem 


+ S 


a 

( 12 ) 


and a similar equation may be written for the p species. In 
Eq. (12) u, V, w are the Instantaneous velocity components In the 
X, y, z directions, respectively, v is the fluid kinematic 
viscosity, and Includes all sources and sinks not directly 

attributable to chemical alteration according to Eqs. (3) and (4). 
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V 


If we assume that an ensemble average of the physical vari- 
ables in the frame of reference is statistically stationary we 
can write any of the variables in terms of its mean and fluctua- 
ting components. If we do this in Eq. (12) and then time average 
the resulting equation according to Reynolds rules ^ we obtain 


ac 




“ + u 


Be 


dx 




Be 


Be 


, 77 Ct 

ST 







e 'u 

a 


B 

“ By 


e 'V ' 

a 



C 


•w ' 

a 


k^(e 


Cft 
a p 


CCA) 

a 


+ s„ (13) 


a 


and similarly for C^ . 

For convenience we assume a statistically steady state (l.e. 
B/Bt =_0) a coordinate_system aligned with the flow direction 
(l.e. V = w == 0 and u = u(z)) , Invoke the two-dimensional 
(B/By = 0) and boundary layer assumptions (B/Bx « B/Bz ) and 
consider no sources or sinks of a (or p) downstream of a contin- 
uous line source of either species. We then obtain the following 
final form of the mass balance equations . 


Be 


— a 


B^C 


a B 


Bz 


« i ^ C'wV - k,(C Co + C'CA) 
2 Bza l^ap ap^ 


(14) 


— 


B% 


p B 


Bz 


^ CAW - ko(C Co + C'CA) 
2 Bzp 2^aP ap'^ 


(15) 


In order to solve Eqs. (l4) and (15)^ we re q uire sev en a d di- 
tiona l parti al diff erential equations for C'w' , Ciw' , C ' 

C'CA j CAT' , CAT' and three additional alg^ralc bquatlofJs for*^ 
a p a ^ P — 2 — p 

the third-order moments C' CA , CACA , C 'CAT . These equations 

are derived In both refs. I and 2 for the to^erested reader. 

In this analysis we have focused on the critical problem of 
the destruction of Co by the NO^ catalytic cycle. Our simulation 
of the NO catalytlc'^cycle Is a simplified one which treats only 
the simultaneous diffusive mixing and reaction of Oo with NO under 
the assumption of the immediate- restoration (and therefore conser- 
vation) of NO by the reaction of NO 2 with 0. The NO^ catalytic 
cycle 

NO + 0 ^ NO2 + Og (16) 


NOg + 0 NO + O 2 


( 17 ) 
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can be approximated for model solution by assuming that NO + NO2 = 
N0;j5- Is conserved and the ratio of NO to NO2 Is a constant In time. 
This Implies the further assumption that the restoration of NO by 
Eq, (17) Is not limited by the supply of 0; a questionable assump- 
tion In the dark sky since. 0 Is produced primarily by photolysis of 
Oo and NO2. This system cannot pretend to follow all of the factors 
which contribute to the quasi -equilibrium concentrations In the 
lower stratosphere. However^ the combined diffusion and simplified 
chemical reaction model used here does capture the main features of 
the NOx catalytic cycle In controlling the local fate of 0^ during 
most turbulent diffusion regimes to be expected In the lower strato- 
sphere. The most questionable assumption is that the production of 
0-^ In the plume Is negligible . 


DESCRIPTION OF COMPUTATION MATRIX 


A total of 213 computations were made using various combina- 
tions of the Input variables to the coupled chemistry/diffusion 
model In order to simulate the reaction of 0^ with NO In the post- 
vortex SST plume. These computations required approximately 500 
hours on the META-4 computer system. The primary Input variables 
which were Investigated were: (l) the turbulence fields (2) the 

source strength and geometry of the Initial species concentrations, 
and (3) the reaction rate constants. In addition, 20 computations 
were made using the three dimensional diffusion model to Investigate 
the effect of the diffusion scale length on plume growth. These 
computations required approximately 100 hours on the META-4. 

Table 1 summarizes the range of variation for each of the Input 
variables. The following paragraphs describe this computation 
matrix . 


TABLE 1 


RANGE OF INPUT VARIABLES 

Variable 


Values 

Turbulence Field 

R1 

0, .1, .2, .35. .5 


50, 100, 200, 500, 1000 meters 

Source Geometry 

mm 

25, 50, 100 meters 

and Strength 

0 

y 

10^, 103, 10^, 3x10^, 105 meters 

Reaction Rate Constant 

^1 ■ 

5x1o5, 5x10^, - 5x1o3, 3x10^ ’ 

Diffusion Scale Length 

^D 

Aiji 

1., 2.5. 5. 10 
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The Turbulence Field 


The turbulence field was simulated by the decoupled Invariant 
turbulence model. The Inputs to this model are the mean wind speed 
and temperature profiles for the altitude Interval of Interest and 
three scale lengths which are derived from best estimates of the 
Integral scale length (Aip) of the turbulent motions. 

The ranges of mean wind shear and temperature lapse rates are 
fairly well known for the_lower stratosphere. For this study, five 
comblnaWons of shear (^u/3z) and potential temperature lapse 
rate (BT/3 z) were chosen to simulate varying degrees of stability. 
The five combinations corresponded to Richardson numbers of 0 
(neutrally stable) to 0.5 (highly stable) and are illustrated in 
Fig. 1. Note that the present model parameters used In the turbu- 
lence model yield a critical Richardson number of 0.55. Therefore, 
all of the flow fields used In this analysis were based on equili- 
brium turbulence- . 

Quantitative Information regarding the Integral scale length 
of the turbulence in the lower stratosphere is meager. In order to 
examine the sensitivity of the chemistry/diffusion model to varia- 
tions In the turbulence field produced by changes In the integral 
scale length, we have chosen five values of the scale length for 
each of the five above-mentioned Richardson numbers. Thus, twenty- 
five (5 Aipx5R-l) different turbulence fields were Investigated. 

Because for the purposes of this analysis, the wind shear and 
temperature lapse rates are functions of Richardson number, the 
equilibrium turbulence field Is completely specified by the combin- 
ation of Ri and Arp -. In Figs. 2 - 5^ we show the sensitivity to 
R1 and Ap of the turbulence parameters which are Inputs to the 
chemlstry/dlf fusion model. 

Q 

Figures 2 and 3 show t he tu rbulent energy (q ) and the verti- 
cal velocity fluctuations (w^) respectively. In each case, there 
Is a strong dependence on both Ri and A,p . Each decreases rapid- 
ly as RI Increases and each Is directly proportional to A^ . 

Figure 4 shows the Reynolds stress (u 'w ' ) . It rapidly ^ 
approaches 0 as RI Increases and is directly proportional to Ap . 

Figure 5 shows the vertical heat flux ( w 'T ' ) . The sensitiv- 
ity to RI Is different from the preceding parameters because the 
heat flux Is 0 not only when RI > Rig = 0.55 tut also when there 
Is no potential temperature gradient (l.e. RI = 0). For a given 
RI , the heat flux Is directly proportional to A^ . 

Table 2 shows the breakdown of the computation matrix for each 
combination of RI and Am . The computations were about equally 
distributed with respect to RI with 42^ of the runs In the highly 
stable regime (RI > .35) and 49^ In the neutral or nearly neutral 
regime (RI < .1) . Approximately 80^ of the runs were for scale 
lengths In the range 100 < < 500 meters . 
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POTENTIAL TEMPERATURE AND VELOCITY 

GRADIENTS 



Richardson number 


Figure 1. The combinations of wind shear and potential 
temperature lapse rate which were used to 
simulate varying degrees of stability In the 
lower stratosphere. 
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TABLE 2 


BREAKDOWN OF RUNS 

FOR Ri AND 

Aiji 


Richardson 
Number (Ri) 

Turbulence Scale Length (Aij.) 

Totals 

50m 

100m 

200m 

500m 

1000m 

0 

2 

12 

2 

7 

2 

25 

.1 

9 

32 

13 

20 

4 

78 

.2 

2 

8 

2 

5 

2 

19 

.35 

5 

30 

4 

7 

9 

55 

.5 

9 

16 

2 

7 

2 

36 

Totals 

27 

98 

23 

46 

19 

213 


Table 3 shows the run breakdown of the turbulence fields in 
terms of decades in Note that the analysis was performed for 

a five decade range of q2 and that about k 0 % of the computations 
were in the interval 10"3 < q^ < 10"^ meter2/sec^. 


TABLE 3 


BREAKDOWN OF RUNS 

FOR q^ 

P 2 2 

q Interval (m /sec ) 

•Number of Runs 

OJ 

1 

0 

1 — 1 

1 

00 

1 

0 

1 — I 

25 

-2 -1 
10 - 10 

62 

10'^ - 10° 

68 

1 — 

0 

i-H 

1 

0 

0 

1 — 1 

43 

1 2 

10-^ - lO'^ 

15 

Total 

213 
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Source Geometry and Strength 


We assumed that the Initial NO concentration distribution was 
blvariant normal and therefore could be specified by a central 
value and horizontal and vertical variances, and q 2 , in the 

following manner ^ ^ 


C^(z,y) 


a., 


.1 


a a 

y z 


exp 



(18) 


where aj^ is a constant to be evaluated from the initial source 
strength' for NO. For the purposes of this analysis, as an initial 
condition the NO was symmetrically distributed about the axis of 
the plume with a standard deviation of 25 meters in both directions 
(l.e. Gy = = 25 m) and a centerline concentration of 2xl0“°ppp. 

The Oo was assumed to be Initially distributed in a complementary 
gaus Sian profile with an ambient value, outside the plume, of 
2x10 °ppp. For these conditions, a^ = 1.25xl0"3ppp m" . 

The present coupled model considers diffusion in the vertical 
direction only. The initial concentration profile is given by 


C(z) 


Cj_(0,0) exp 


1 

z 


(19) 


and Cj_(0,0) is obtained from Eq. (l8) for each combination of a-y 
and Gg . ^ Table 4 summarizes the values of the initial centerline‘s 
concentrations for each of the source geometries , 


TABLE 4 


INITIAL CENTERLINE CONCENTRATION (ppp) 

G 

Z 

(meters ) 

G (meters) 

Vi 1 ^ 

OJ 

o 

1 — I 

io3 

o 
1 — 1 

3x10^ 

105 

25 

50 

100 

.5x10“^ 

.25x10“^ 

. 125 x 10 "^ 

.5x10-'^' 

... 25 x 10 "'^ 

. 125 x 10 "^ 

.5x10-3 

.25x10-3 

.375x10‘9 

. 15 x 10-3 

.75x10"^ 

.375x10-9 

.5x10“9 

.25x10-9 

.125x10”9 


Thus, a combination of Oy and a may be related to a 
particular time in the diffusion of a plume whose source strength 
is fixed. 
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Alternatively, because only enters the coupled model 

directly, a„ may be used as a measure of the source strength as 
follows. Consider the total amount (Q) of NO In a vertical 
column per unit distance of travel and per unit width of the plume. 


Q = J C(z ) dz 
— 00 

If we Integrate Eq. (20) using Eq. (19)^ then 


( 20 ) 


But, by Eq. (l8) 


and, therefore. 


Q -- v/2^ a C . (0, 0 ) 
z ,x 


a. 


C. (0,0) = — — 
' o o 

y z 


Q = 


a 

y 


(21) 


( 22 ) 


(23) 


Because Is constant, Eq. (23) demonstrates that the amount of 

NO In a vertical column (l.e. the source strength) Is Inversely 
proportional to Cy . Thus, If we compare computations at a fixed 
time after emission, then cy determines the source strength and 
determines the distribution of NO for a fixed strength. 

In Table 5, we summarize the number of runs used for each 
combination of a and . Approximately one half of the runs 
were for Oy ■= lOO meters. Less than 5^ of the runs were for 
a > 10^ meters. 

■ u 

TABLE 5 


BREAKDOWN OP RUNS FOR 0 ^ and 0^ 


Vertical Spread 

Horizontal Spread 

~ a (meters.) 

^ L 

_ Total- 

a (meters) 

OJ 

0 

1 — 1 

103 

10 ^ 

S xld^ 

105 


25 

55 

12 

21 

0 

0 

88 

50 

28 

14 

17 

•4 

6 

69 

100 

31 

6 

19 

0 

0 . 

56 

Total 

114 

32 

57 

4 

6 

213 
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Approximately 4o^ of the runs were for = 25 m and the remain- 
der of the rims were about evenly divided between :o^ = 50 m and 
= 100 m . 


Reaction Rate Constant 

We have chosen a simplified version of the NO^ catalytic cycle 
for the destruction of Oo. We have assumed that the rate of 
production of 0^, In the exhaust plume during the early post-vortex 
stage is negligibly small compared to the destruction of Oo by NO 
and that photolysis of NOg and the abundance of 0 In the plume are 
sufficient to restore all NO 2 to NO + Og. Thus, the catalytic 
cycle is simply 

NO + 0^ NOg + Og (24) ■ 

NOg + 0 NO + Og ( 25 ) 


and 


a 



chem 


k^Qo^llNO] + [0^] '[NO] '] (26) 


a 


[NO] 


chem 


0 


(27) 


The sensitivity of the model to kn was Investigated by vary- 
ing k^^ through four orders -of -magnitude from 5xlo5 to 5x10^ 

(pp - sec)"4__ Table 6 summarizes the number of runs for each k-j^ . 


TABLE 6 

BREAKDOWN 

OF RUNS FOR k^ 

k f ) 

1 pp-sec ' 

Number of Runs 

5x10^ 

122 

■ 5x10^ 

31 

5x1o3 

31 

5x1o2 

29 

Total 

213 
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About 60 % of the runs were with k]_ = 5x10^ 1/pp-sec. and the 
balance of the runs were equally distributed among the other three 
values of k-]_ . 

Diffusion Scale Length 

Extensive experimentation has shown that during the initial 
phase of plume diffusion the diffusion scale length is correctly 
determined by the width of the diffusing plume. At some point, 
however, this scale length becomes constant and the diffusion 
process goes over to the classical Plckian rate. It can be shown 
that the upper limiting value of the scale length (Aj^) is propor- 
tional to the turbulence scale (Aj,) although the constant of 
proportionality is poorly defined. ¥e have examined four ratios of 
Aj)/Ap (l^ 2.5^ 5 and 10) for each of five turbulence fields. 

Since the NO is conserved by the NO^ catalytic cycle, the 
diffusion of NO can be decoupled from the chemistry. Therefore, 
the sensitivity of the plume growth to changes in Ajy^Am could be 
investigated using our three-dimensional diffusion model. The 
results of these 20 computations together with the results of the 
213 computations using the coupled chemistry/diffusion model are 
discussed in the following sections. 


DISCUSSION OF SENSITIVITY ANALYSIS 


In. this section, we present the results of a sensitivity 
analysis of the two-dimensional coupled chemistry/diffusion model 
(as well as the three-dimensional diffusion model) to the primary 
Independent variables of the model. We will concentrate on the 
principal output parameters of the models such as plume growth, 
total reaction rates, etc. Keep in mind, however, that by emphasiz- 
ing these results, we will, of necessity, omit a wealth of data 
relating to the details of the computations. For example, the time 
histories of the profiles of all turbulent correlations and the 
various balances which exist between the fluxes of these correla- 
tions in the dynamic equations are available but are not shown. 

These details can be described at a later date when more detailed 
measurements of jet wakes become available. 

Plume Growth 

In our model of the NO^ catalytic cycle, the NO is conserved 
and therefore acts as a surrogate tracer for any nonreactlve exhaust 
product. Hence, we can separate the effects of turbulent diffusion 
and chemistry by following the growth of the NO plume. To do this, 
we have calculated the three -dimensional diffusion history of a 
plume, for which the initial horizontal and vertical variances are 
625 m2, over a time Interval corresponding to a concentration 
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Vi 


dilution of nearly 10^. These calculations were made for five 
equilibrium turbulence fields for which the Intensity of turbulence^ 

, "varied from 0.008 to 0.5 m^/sec^^ r 1 from 0 to 0.5^ and the 
scale length of the turbulence^ Aiji , was 100 meters. 

The Independent variables of the diffusion model are the mean 
wind shear and potential temperature lapse rates (which together 
define R1 .) , the Intensity of turbulence^ the vertical and lateral 
velocity fluctuations, and the vertical heat transfer. The sensi- 
tivity of each of these variables to R1 and to Aip was previously 
described In Pigs. 1-3 and 5. 

Besides these Inputs, the diffusion model requires specifica- 
tion of the diffusive scale length, Ap . This latter length Is 
analogous to, but not necessarily the same as, the turbulence scale 
length because Ap depends on the dimensions of the plume. Prior 
applications of this model have shown that the diffusive scale 
length for Individual plumes Is very well specified by the standard 
deviation of the concentration distribution at least until this 
length reaches a multiple of the turbulent scale length. This 
multiple, however. Is not well defined. We have chosen four values 
of Ap/Aip , ranging from one to ten. In which to Investigate the 
sensitivity of plume growth ( Cy and a|) to diffusive scale length. 
Keep In mind that both and"^ are Inputs to the coupled model 

and that, therefore, the following discussion will also provide a 
frame of reference when we get to the coupled model results. 

Figure 6 Illustrates a typical case of the diffusion history 
of a plume. It shows the lateral and vertical variances of the, 
calculated distribution, and , and the ratio of the concen- 
tration at the center of the plume to Its Initial value as a func- 
tion of the t,lme of travel. In a time Interval of lo5 seconds 
(approximately 28 hours), the plume width, as measured by Cy and 
has spread approximately two orders of magnitude In both the lateral 
and vertical d,lrectlons, and the concentration has been diluted by 
a factor of 10^. Because the lower stratosphere Is stably strati- 
fied, the lateral spread Is larger than the vertical spread. The 
degree of anisotropy as measured by ^ = 1^)5 is 0.2 . 

Although the data are not, shown, tbe anisotropy Increases (l.e., 
decreases) with Increasing R1 or decreasing q'^ . In fact, for 
q2 < 10"2^ the diffusion. Is almost entirely In the lateral direction. 

The change In both and proceeds more rapidly than f*"! 
until either c2 or (or both) reaches the value of the cut.off 

Aj^ , which In this case Is = 6.2x10^ m^ . After this point Is 
reached, the plume spread goes over to t+1 In the lateral direc- 
tion (a result expected from Plcklan diffusion theory) and to less 
than t+l In the vertical direction. 

We could present similar plots for each of the other 19 cases 
which were Investigated, However, we have found that all results 
could be conveniently summarized In one nomogram. This graphical 
tool not only summarizes the results but also yields the partial 
sensitivity of a computed parameter to each of Its Independent 
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variables. Because we will use the properties of this nomogram 
frequently, we will digress here and discuss Its properties. 

Figure 7 Illustrates the generalized nomogram. Consider a 
dependent variable P to be a function of n dependent variables 

Xn, Xg, Xj^ (where n = 4 In Fig. 7). One may easily determine 

the value of P for any combination of Xj^ by starting In the 
upper left quadrant with Xn and marching around the diagram to 
P . For example, the variables x^ and Xg define a point on the 
X2) axis, fo and Xo define a poxnt on the fo(x]^, Xg^ 
Xo) axis, and finally fo and ^4 define a value of P . If 
there were additional Independent variables X(- , xg , etc., one 
would simply continue around the diagram untll^the P axis was 
reached. 

When the scales for each of the axes of the nomogram are 
logarithmic, the nomogram yields the partial sensitivity of P to 
any fj^ . To demonstrate this property, assume that ^ 

^2 } ^ 2 ) ^ logarithmic scales In Fig. 7. Now consider the 

Increment In P produced by an Increment In x^ as shown by the 
two dashed curves. If we use the definition for the slope of a 
line, then 

A in fg 
A in x^ "" ^1 

where S;^ 4s the slope of the xg = constant curve. Note that we 
need not limit the analysis to straight lines If we utilize the 
slope of the tangent to an Iso-xg contour. If we continue around 
the diagram, then 


and 


therefore , 


or 


A in f^ ^ 
A in fg ^ 


A in P 
A in f^ 

A in P 
A in x^ 

A in P 
A in x^ 


/ A in fg \ / A in fg \ / ^ P 
yA in x^y y A in fg y y A in f^ 


In the limit of small changes In x^ , we have 

a in P _ ^1 ^3 
d in x^ Sg 


(28) 


(29) 


where , the slope of the tangent to the curve of constant 

Xj_^]^ , Is an amplification factor and d in P/d in x^ , Is the 
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partial sensitivity of P to . A similar analysis would show 
that the sensitivity of P to fg Is simply 

h £n P _ fs 
d Jin fg Sg 

and so on for any function fj_ . The final sensitivity of P to 
Is not given directly but can be obtained by noting the Incre- 
mental change In P for a change In while, holding f~ 

constant. 

The result given by Eq. (29) can be generalized to the case of 
n Independent variables as follows 


^ in P . . . . 

a in f^. - Sg Sg . . . , 


(30) 


where = 1 for k > n and Is any transfer function 

(fj, = f]_j fg^ or f^ In our Illustration). 

We can see Immediately from Eq. (29) that If either Sn or-.- 
So Is 0 (l.e. contours of xg or x^ are horizontal) or 89 Is 
00-^ (l.e. contours of xo are vertical) , then x. will have no 
effect on P . If all the are nonzero^ then the largest value 

of Is the dominant Influence on P . Finally, the algebraic 

sign of Sj^So/Sg determines whether P will Increase or decrease 
with a glven'^lncremental change In . 

Nomograms which summarize the diffusion computations are shown 
In Figs, 8 and 9 for Cy and respectively. The Independent 

variables which coincide with x^, Xg, x^ and x^^ are R 1 , , t , 

and . We can, of course, recover the results previously shown 

In Fig. 6. For example. If we choose R 1 = 0 . 1 ^ and A^ = 100 m, 
then q 2 = fg = 0.24 mvsec^ , and when t = 10^ sec ana Aj) = 250 m 
we get Oy = 680 m . 

In Fig. 8, the a curves coincide for o < Ap . Beyond the 

cutoff point, the plume spreads more rapidly for larger Ap . 

However, the effect of Aj^ on a Is significant only for rela- 
tively large turbulence Intensities and for very large times. For 
example, when q^ = 0.1 m^/sec^ and t = 10^ sec , a = 3OOO m 

when Aj) = IQOO m , but Cy Is only 1000 m If Aj) = 100 m . For 

a smaller q^ = 0.001 , the difference In Is 300 m relative to 

250 m for the same range of A^, and at the same time. For the 
latter Intensity of turbulence, noticeable differences In plume 
growth due to Aj^ effects would only occur In a time scale of 
weeks . 

The sensitivity of Oy to A^, Is summarized In Fig. 10 where 
contours of constant f^ are shown currespondlng to Fig, 8. Only 
when f^ Is large does"^ A^, influence a . When f^ Is small. 
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The dependence of the horizontal spread of the 
concentration distribution, Cy , on the turbu- 
lence field (R1 and A>p) and the diffusive scale 
(Aj)) , for ticaes after emission between 15 min. 
and 28 hours . 





Figure 9. The dependence of the vertical spread of the 

concentration distribution, , on the turbu- 
lence fields (R1 and Ajp) and the diffusive scale 
length (Aq) , for times after emission between 
15 min, and 28 hours. 






Ap(m) 

The sensitivity of the horizontal sprea^ 
to the diffusive scale length, Aj) , fo: 
values of fo where fo is a function 
Aij. and t , 




changes, in Aj^ do not affect the plume spread. As noted in Fig. 8^ 
the magnitude of , f ^ increases with t and q^. 

Because of vertical stability^ can be considerably smaller 

than a . As noted in Pig. 9, the effect of on is much 

smaller^ than was the case for Oy . 

Figures 8 and 9 demonstrate quite clearly the sensitivity of 
both Cy and to the turbulence field. In particular, they 
demonstrate the sensitivity of Cy(and o^) to R1 . We could 
have just as easily used = Ap in our nomogram with the same 
results. The amplification factors (Sp, 82^ Sq in Eq. (29)) . show 
that the turbulence field is the critical parameter in determining 
the growth of the plume . For the entire domain of the diagram. 

So = 1 and So < 0.5 . However, Sp approaches -00 (the minus 
sign denotes a'^decrease in a for an Increase in Rl) as the 
stability Increases and approaches the critical value. At Rl^ , 
of course, the diagram is no longer applicable since equilibrium 
t.urbulence cannot be sustained. Because experimental evidence does 
suggest high stability in the lower stratosphere. Fig. 8 (and 9) 
demonstrates that one had better know the turbulence before he 
attempts to compute the dilution time history of the plume . 


In Fig, 11, we further demonstrate the Importance of the 
turbulence intensity by showing the time required to reach a concen- 
tration dilution of 103, t^ . We postulate that such a dilution 
would be the limit of detectability of the plume and therefore t^ 
is a good first estimate of the duration of Phase III. It is clear 
that t^ is much more dependent on q^ than on Ap, . ■ An order of 
magnitude change in q^ also changes t^ by an order of magnitude, 


whereas an order of magnitude change in Ap, changes t^ 


on the 


scale of hours. When q^ < 10~2 m^/sec‘^ measured in weeks 

and therefore effects are insignificant. When q2 > o.l, Ap, 

affects are larger but the results are really of academic Interest 
since the plume will be completely diffused into the environment 
within hours anyway and furthermore the probability of sustained 
turbulence of that intensity is slight. 


Chemical Reactions Between NO and 0^ 

In this analysis, we concentrate on the destruction by the NO^ 
catalytic cycle of only that component of Oo which occurs in the 
ambient (nonplume) environment. It is recognized that the NO in an 
Intact plume also reacts with and destroys locally produced Oo. 

This local production of Oo may be high because of excessive produc- 
tion of 0. However, the NO concentrations are also high and the 
locally produced Oo is destroyed in a "do-nothing" cycle. In ref 2, 
it was shown that bhls cycle persists until plume dilutions of lo3 
are achieved. It can therefore be assumed that there is a limit- 
ing value of the NO concentration just capable of destroying 0^ 
at its local production rate and any excess of NO above tMs level 
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EFFECT OF TURBULENCE INTENSITY 
ON DURATION OF WAKE PHASE 
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2 

Figure 11. The effect of the turbulence intensity, q , and the diffusive scale length 
Aj) , on the time required to achieve a concentration dilution of 103 , t^ . 
This time is an estimate of the duration of the SST wake phase. 




is inactive. 

The depletion rate of ambient 0^ depends on the relative rates 
of diffusive mixing and chemical reactions , For a vertical column 
of unit cross-section and extending through the plume^ the actual 
depletion rate of ©^(R) Is given by 

/ oo 

— 00 ^ 


where a denotes the Oo and P denotes the NO, The maximum value 
which C can attain In the plume is the environmental concentra- 
tion of Oo and In this case 0 ^ is zero because the 0 ^ Is uniform- 
ly distributed. The maximum depletion rate of 02 (R*) -ls therefore 
given by 

/ OO 

O-p dz (32) 

e -00 


where Co.^ Is the ambient concentration of O 3 and the Integral is a 
measure of the total amount of NO in a vertical column through the 
plume . The chemical efficiency for any combination of NO and 0^ is 
then given by the ratio of Eqs. (31) and (32) 


E 


R* 





C'CA) dz 
g p' 



dz 


(33) 


There are two effects which prevent the reaction from proceed- 
ing at Its maximum rate. Firsts the ambient Oo may not be diffused 
Into the plume rapidly enough. That Is, the turbulence Intensity 
may be too low to diffuse Oo into the plume after the ambient Oo in 
the plume Is destroyed. Because there Is more than enough NO to 
react with locally produced O 3 , this means that substantial amounts 
of NO may be inactive until plume dilutions of at least 103 are 
achieved. It was shown In a previous section that dilution of this 
magnitude may not occur for days to weeks for typical turbulence 
intensities in the stratosphere. The diffusion limitation may be 
quantitatively Investigated by defining the parameter, Lj^ , as 
follows p°° 





dz 


C Co dz 
a p 



(34) 


In this form, Lj) represents the deficit in the reaction rate due 
to slow diffusion of ambient 0^ . 
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The second limi tatio n to the maximum depletion rate is due to 
the mlxedn£ss term . Neglect of this term is equivalent to 

assuming and Co are constant In time - a situation which 

certainly does not describe the real situation at least during the 
early stages of the reaction. Quantitatively, the mixedness limit- 
ation, Ljy[ , Is defined as 





CCA dz 
g p 



(35) 


It should be noted that Eqs . (33) through (35) are simply tied 
together by the relation 


E = 1 



(36) 


The Invariant coupled diffusion/chemistry model has been 
exercised In order to demonstrate the sensitivity of E (really 
Ljj and L]v[) to the various Independent variables for a wide range 
or equilibrium turbulence fields, plume geometries, source strengths 
and reaction rate constants. Because the model In Its present form 
can only handle two-dimensional diffusion. It Is Implicitly assumed 
in the following discussion that the lateral extent of the plume is 
large relative to the vertical extent and that lateral concentra- 
tion gradients are negligible compared with the vertical. The 
three dimensional diffusion calculations which were described 
earlier demonstrate the validity of this assumption. 

All of the computations are summarized In nomograms which will 
be discussed shortly. Because of the complexity of the computa- 
tions and the number of Independent variables which are Involved, 
the Information contained In the nomogram may not be Immediately 
obvious to those who are more familiar with conventional sensitiv- 
ity plots where one Independent parameter is varied and all other 
parameters are held constant. It is therefore beneficial to 
describe a few of the Individual computations In order to under- 
stand the underlying physics which controls the Interaction between 
turbulent diffusion and chemical reactions and in order to s.et the 
stage for the nomograms. 

Keep in mind that all coupled model computations are dynamic 
and that time (or more correctly for a shear flow the distance 
downstream of a steady line source) Is a variable and the results 
must be qualified by stating the time regime over which the results 
are applicable. Figure 12 shows a typical time history of the 
total reaction rate as a function of the turbulence intensity. 
Initially, the reaction proceeds rapidly as the Oo present In the 
plume reacts with the NO. While the local 0^ Is being depleted, 
the reaction rate steadily decreases until It just balances the 
diffusion rate. At this time and for some time afterward, a quasl- 
equllibrlum period exists during which turbulent diffusion provides 
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REACTION RATE TIME HISTORY 
Richardson No. =0.2 
CTv(0)=l00m CT7{0) = 25m 
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Figure 12. Typical reaction rate time histories for five values of turbulence 

Intensity varying from 0.03 to 12.8 m2/sec^. R represents the integral 
over the vertical extent of the plume of the instantaneous local reaction 
rates at time t . 



0^ as rapidly as the chemical reaction depletes It. After the 
equilibrium period^ the reaction rate drops as the result of 
further NO dilution. 

Each computation was taken at least as far as this equilibrium 
period. The results presented below are for a distance of 50 km 
(approximately 15-30 minutes) downwind of a steady source whose 
Initial concentration Is given In Table 4 for various combinations 
of Cy and . Recall that the Initial concentration for the 
coupled model Is Inversely proportional to the product c7yO (see 
Eq. (l8)). Thus, the equlllbrl\im period can be related to Phase III 
time using the three-dimensional diffusion nomograms described 
earlier. In the following paragraphs, the separate effects of turbu 
lence, source geometry, source strength, and reaction rate constant 
on the species concentrations and reaction rates are described. 
Implicit In these results Is the assumption of quasi -equilibrium. 


The Turbulence Field 

Consider, first, the effect of varying the turbulence field on 
the NO-O^ reaction. In this discussion, the Independent variables 
are R1 "^and Am . The turbulence Intensity, q^ , Is In effect a 
dependent variable and Is therefore not shown. We could just as 
easily have used R1 and q^ as the dependent variables without 
changing the basic results but the resulting curves would have 
masked the sensitivity of the model to the turbulence scale length 
The three variables are, of course, simply related for equilibrium 
turbulence as was shown In Fig. 2. 

Figure 13 shows the effect of R1 and A/p on the ratio of 
the maximum concentration of NO to Its Initial value. The reduc- 
tion of NO at the center of the plume Is solely the result of 
diffusion because the NO Is conserved and Is unaffected by the 
chemistry. It Is evident that diffusion proceeds more rapidly as 
Arp Increases and as R1 decreases (l.e., as q^ increases). 

Figure l4 shows the effect of the same turbulence fields on the 
reaction rate. A two-order of magnitude Increase In the reaction 
rate can be observed as the turbulence field varies from highly 
stable, small scale eddies to neutrally stable, large scale eddies. 

Figure 15 shows the 0^ concentration at the center of the 
plume. It Is evident that when the turbulent diffusion rate Is low 
the ©2 concentration can be more than three orders of magnitude 
lower'^than Its ambient value. In other words, a portion of the 
plume Is essentially free of Oo - a region which Is termed the 
"Oo - hole.". The presence of the Oo hole has a detrimental effect 
on'^the chemical efficiency because all of the NO present In the hole 
(and keep In mind the hole Is at the center of the plume where the 
NO concentration Is a maximum) Is unreactlve. 

The extent of the 0^ hole can be quantitatively estimated If 
one Is willing to arbitrarily define a concentration level below 
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[NO]„ax/[NO] 


EFFECT OF TURBULENCE FIELD 
ON MAXIMUM NO CONCENTRATION 



Figure 13. Typical effect of the turbulence field, (as defined 
by Ri and Ajp) , on the ratio of the NO concentration 
at the center of the plume to its initial value. The 
initial plume geometry ( 02 ( 0 ) = 25 m and ay(0) = 
100 m) . is fixed. 
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Normalized reaction rate , R/R 



;ure l4. Typical ' 
depletlo: 
geometry 
rate con 












which the Oo concentration Is effectively zero. For purposes of 
discussion. If we define this level as 10“^ times the ambient 
concentration, then Fig. l6 shows the maximum thickness of the Oo 
hole which Is developed. For very stable and small-scale turbu-'^ 
lence, the thickness of the hole can be more than 100 meters or 
nearly equal to the Initial thickness (l.e,, of the plume. 


Plume Geometry 

The effect of plume geometry was Investigated by varying both 
0y and . Because only vertical diffusion Is considered, the 

lateral extent of the plume Is, In effect, a measure of the source 
strength as was noted earlier In Eq. (23). We will, therefore, 
postpone the a„ discussion until the next section and consider 
only In this section. 

Figure 17 shows the effect of on the maximum NO concen- 

tration for five turbulence fleMs . For a fixed source strength, 
the concentration gradient (^C/Sz) decreases as Increases 

and, therefore there Is less diffusion of NO from the center of the 
plume for the larger values of . 

Figure 18 shows the Oo concentrations at the center of the 
plume as a function of for three turbulence fields. For a 

neutrally stable environment, the Oo concentration slowly decreases 
as 02 Increases. As the stability Increases, the balance between 
the reaction rate and the diffusion rate shifts and the sensitivity 
of Oo concentration to 0£ Increases. Although the sensitivity Is 
Increasing with R1 , the concentration levels are 10“3 to 10~5 
times the ambient value and therefore we can anticipate that these 
results will have a minor Influence on the chemical efficiency. 

Figure 19 shows the effect of O2 on the total reaction rate 
for the same three turbulence fields as Fig. 18. For the most part, 
the reaction rate Is insensitive to o^. and never varies by more 
than 20^ for a factor of four change In 0 . Again, these results 

suggest that 0^ will have a negligible effect on the chemical 
efficiency. 


Source Strength 

Figure 20 shows that the normalized maximum concentration of 
NO Is Independent of source strength or. In other words, the 
maximum local concentration of NO varies directly as the amount of 
NO emitted at the source. This is not surprising because the NO 
Is conserved and the Initial concentration Is directly proportional 
to the source strength. 

The source strength has a significant effect on the Oo concen- 
tration at the center of the plume as shown In Fig. 21. For a fixed 
turbulence field, as the amount of NO Is Increased (l.e. Cy 
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EFFECT OF TURBULENCE FIELD ON 
THICKNESS OF O3 HOLE 

k'SxlO® (l/pp-sec) 
o"y (0) = IOOm 



Figure 16. The effect of the turbulence field on the 
maximum vertical thickness of the Oo-hole 
for the same conditions as for Fig. l4. The 
Oo-hole Is defined as the region around the 
center of the plume within which the 0^ 
concentration Is less than 10"3 times the 
ambient 0^ concentration of 2 xl 0 “O ppp. 
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;ure 17. Typical effect of o^(0) - the Initial vertical 
spread of the NO concentration distribution - on 
^the ratio of the NO concentration at the center 
'of the plume to its initial value for five yalues 
of R1 . The parameters /bp , Oy(0) , and x ar 
held fixed at 100 ra , 100 m , and 50 km respect Iv 








Normalized reaction rote, R/R 








EFFECT OF o-y(O) ON 
MAXIMUM NO CONCENTRATION 


cTz (0) = 25m 



o-y (0), m 

Figure 20, Typical effect of Cy(0) - the initial lateral spread of 
the NO concentration distribution - on the ratio of the 
NO concentration at the center of the plume to its init- 
ial value for five values of Ay . The parameters Rl- , 
o^{0) , and x are held constant at 0.1^ 25 m , and 50 km 
respectively. For two-dimensional diffusion, Cy(0) is 
Inversely proportional to the source strength. 
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EFFECT OF o-y(O) ON 
MINIMUM O3 CONCENTRATION 


k=5xlO^ (l/pp-sec) 

CT2(0)=25m 



o-y(O), m 

Figure 21, Typical effect of ay(0) on the Oo concentration 

at the center of the"^ plume for the same conditions 
as for Rig. 20 and for k = 5xl05(p.p-sec ) . 


decreases) the total reaction rate increases and therefore the 0 ^ 
concentration decreases. The sensitivity of [03]fjj^n 
Increases as the turbulence Intensity (i.e., Ajj,) decreases 
because the O 3 is less able to penetrate deep into the plume. That 
is to say, the Increased sensitivity Illustrates the tendency for 
the reaction to become diffusion limited for low turbulence Inten-^ 
slty. 

Figure 22 shows the effect of source strength on the reaction 
rate. For very dilute plumes, the reaction rate should be linearly 
dependent on the amount of NO present in the plume. The curves 
shown in Fig. 22 appear to be approaching this linear dependence 
on a for a > 10^. As the amount of NO is Increased 
(dy <^ 10 ^), the reaction rate departs considerably from a linear 
dependence due to diffusion and mixedness limitations. That is 
to say, diffusion and mixedness are becoming increasingly important 
in controlling the total reaction rate. For example, a tenfold 
Increase in the amount of NO (from Cy = 103 to o = 10‘2) in a turbu- 
lence field whose = .06 m^/sec^ \Aqi - 50 m) results in a 
negligible Increase in the reaction rate. This is because the 
diffusion rate is extremely- slow, the reaction is occurring at the 
edges of the plume and the bulk of the NO has no 0- with which to 
react. ^ 


Reaction Rate Constant 

The reaction rate constant for the depletion of Oo was varied 
through four orders -of -magnl tude . Because the NO is conserved, k 
has no effect on the NO concentration. It does have, however, a 
pronounced effect on the O 3 concentration at the center of the 
plume as shown in Fig. 23 and on the reaction rate as shown in 
Fig . 24 . 

In each figure, for large values of k both the minimum Oo 
concentration and the reaction rate are strongly dependent on the 
turbulence intensity. That is to say, the reaction is tending 
toward the diffusion limitation. For small values of k , the 
reaction is nearly independent of the turbulence intensity and 
proceeds at a rate controlled by k , The reaction is now slow 
enough so that turbulent diffusion can supply all of the O 3 needed 
to react with the NO and, in fact, the Oo concentration doesn't 
vary by more than a factor of two across'^the plume. Note in 
Fig. 23, that a tenfold decrease in k from 5xl05 to 5 x 1 Q 4 
practically eliminates the Oo-hole. Also note, this time in 
Fig. 24, that the total reaction rate increases as k^ and that 
n = +1 for small" k and that n decreases steadily with increas- 
ing k . We would expect, therefore, that the chemical efficiency 
would decrease with Increasing k because the maximum possible 
reaction rate is everywhere proportional to k^ . 
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Normalized reaction rate, R/R 


EFFECT OF cry (0) ON 
REACTION RATE 

(Tz (0) = 25 m 
Ri = O.I 
x = 50 Km 


1000 m 


R*=2xl0 


Figure 22. 


Typical effect of a (0) on the total depletion 
rate of ambient 0^ for the same conditions as for 
Fig. 21. 



EFFECT OF k ON 
MINIMUM O3 CONCENTRATION 


Ax=IOOm 
CTy (0)=I00m 
CTz (0) = 25 m 



\ 0 ~^ 10 "' 10 ° 


k/K 


Typical effect of the reaction rate constant on 
the Oo concentration at the center of the plume for 
five values of R1 . The other parameters are held 
constant; TW, = 100 ra , 0 ^( 0 ) = 23 m , o (O) = 100 
and ' X = 50 m . 
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Chemical Efficiency 


The results described above can be recast in terms of the 
chemical efficiency (E) and in terms of the diffusion and mixedness 
limitations (Lj) and L]y[, respectively) which act independently to 
limit E to values less than one. In an earlier section, each of 
these terms was defined (Eqs. (33)-(35)). 

One might well recall, at this point, the implications of the 
chemical efficiency as far as the Oo inventory in the stratosphere 
is concerned. In Ref, 2, we showed'^that the average chemical 
efficiency for the entire residence time of the plume E (Tj^) is 
given by 

E(T^) = 1 - f- (1 - E) (37) 

R 

where T is the duration of Phase III and Tp is the , residence 
jtlme of the plume in the stratosphere. Clearly, as T/Tp -*• 1 , 

E (Tj^) -► E and E , therefore, is a measure of the overall effect 
of Phase III chemistry limitations on the theoretical depletion 
rate of environmental Oo due to the SST, Of course, if T/T^ — *■ 0 
or E — ► 1 , then E (Tj^j -► 1 and Phase III chemistry effects are 
of academic interest. Present best estimates of T/Td vary from 
less than 0.1 to ’.as much as 0.5 and as will be shown below E can 
vary from less than 0.01 to 0.99. Thus, Phase III chemistry limit- 
ations can have a significant influence on the Oo inventory in the 
stratosphere. 

The chemical efficiency nomogram is presented in Fig. 25. The 
independent variables are Rl, Ap, Oy and h . All of the data in 
Fig. 25 are for a„ = 50 meters. Before we discuss these results, 
we should note that E is relatively insensitive to . 

The effect of on E is shown in Fig. 26. A band is shown 

for each of three values of Oy . The upper value of the band 
usually corresponds to = 100 m and the lower value to Cz = 25 m. 
For dilute plumes, is seen to be a relatively small effect in 

comparison to changes in Oy . For compact plumes, the spread is 
larger but the differences are between a fraction of and 3^ - 
differences which are of little consequence. We, therefore, omit 
a from further discussion and note that, strictly speaking, all 
or the following results are for = ^0 m . 

Let us return to Fig. 25 in which some of the results described 
earlier are clearly evident. When the turbulence intensity is high 
and a is large, E approaches unity no matter what the reaction 
rate constant. In contrast, when q^ is low and the plume is 
compact k is an extremely important parameter and the efficiency 
can vary from 80^ for slow reactions to less than 1% for very fast 
reactions. 

In the context of our earlier discussion of the nomogram, the 
slopes of the curves in Fig. 25 amplify the sensitivity of E to 
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EFFECT OF CTzlO) ON 
CHEMICAL EFFICIENCY 


<r„(0) = IG'’m 


(0)=l0^m 


<r„(0)=IO^m 
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Figure 26. 


The effect of o^{ 0 ) - the initial vertical 
spread of the NO concentration distribution 
on the chemical efficiency for three values 

of o (0). 


1 







\ 


the R1 . The slope of the Iso-A/p contours (S^^) Increases rapidly 
as the Rl approaches the critical R1 of 0,55, Thus E is 
most sensitive to the turbulence field when the stratosphere is 
highly stable. The slope of the iso-Oy contours is nearly constant 
for less than 0,1 and therefore the sensitivity of E to the 

turbulence field is relatively independent of plume geometry. For 
very large q2 , S2 approaches infinity and E becomes independ- 
ent of the turbulence. However, this result is mainly of academic 
Interest because sustained turbulence in excess of 10 m2/sec2 in 
the stratosphere is highly unlikely. The slope of the Iso-k contours 
(So) is maximum for Intermediate values of the transfer function 
fo'^. For large values of fo (l,e. very compact plumes). So 
approaches 0 and E becomes'^ Independent of Rl and varies only 
with k . Thus, E is most sensitive to the turbulence field when 
the flow field is highly stable, the reaction rate is fast and the 
plume is moderately diffused. 

The partial sensitivity of E to Cy is illustrated in Pigs. 
27-29 for three values of the turbulence intensity and. three values 
of k . The three figures clearly show that: 

1) E Increases strongly with Cy or, in other words, the 
less NO emitted per unit distance of travel the more efficient the 
destruction of 0^ . 

2) E is most sensitive to Cy when the NO concentrations are 

high. As the plume becomes diffused, the sensitivity to Cy 
decreases. ^ 

3) The sensitivity of E to Cy decreases as k decreases. 
Figure 29 shows that E is always above 4o^ when k = 5xlo3 
(l/pp-sec). In this case, conventional mean-value chemistry would 
predict the depletion of 0^ to within a factor of two. 

The sensitivity of E to k is Illustrated in Fig. 30 in 
which lines of constant f^ are shown. Recall that fo is a 
transfer function of Rl, '^Aip and a and that for Increases in 
each of these variables fo Increases with Ri and decreases with 
Ajji and Oy . When the reaction rate constant is very small, E 
is nearly mlty no matter what the value, of Rl , Am and o„ , As 
k Increases, the reaction becomes more diffusion and mixedness 
limited and the efficiency decreases. For very large values of f^^ 
the Iso-f^ contours tend to merge and E becomes a function of k'^ 
only , 

In order to put all of the above results in perspective, we can 
illustrate the Importance of each of the Independent variables by 
citing a typical example. For purposes of discussion, assume that a 
typical condition in the lower stratosphere is given by q2 = 
sec2 ^ Oy = lo3m , and k = 5x10^ (l/pp-sec) . For these conditions. 
Fig. 25 snows that E = 10^ and we can ask the question, how much 
does E change if each of the variables is changed Independently 
by one order of magnitude? Table 7 summarizes the answer to this 
question. 
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Figure 29. The sensitivity of the chemical efficiency to cTy(O) for three 
values of the turbulence intensity between 1 and^l0"2 m^/sec^ 
and for k = 5xl03(pp-sec )“4 . 







TABLE 7 



ILLUSTRATION OF SENSITIVITY OF E . 

TO INDEPENfelfc.™ '■ . : 

(BASE VALUE 0F‘ e' = 10^) '"'-i 


Independent 

Direction of Order 

Value of Variable 

E(%) 

Variable 

of Magnitude Change 



p 

Increase 

10'^(m^/sec^> 

45 


decrease 

10'^(m^/sec^) 

2 


Increase 

lO^(ra) 

4o 

^y 

decrease 

lO^(m) 

3 


increase 

5xl0^(l/pp-sec) 

2 

k 

decrease 

5xlo3(l/pp-sec) 

60 


Table 7 shows that a tenfold change in any one of the independ- 
ent variables can result in a five-to sixfold change in the chemical 
efficiency. It can be demonstrated using Eq. (37) that this varia- 
tion in E is equivalent to E values of between 0.5 and 0.8 
when the duration of the wake phase is one half of the plume 
residence time in the stratosphere. In other words, present esti- 
mates of the impact of the SST on the Oo Inventory must be reduced 
by 20 % to 50 ^ - a rather wide variatlon'^to say the least - for the 
conditions chosen. More accurate predictions will require consider- 
ably better than order of magnitude estimates of the independent 
variables. These results, as was noted above, are for purely hypo- 
thetical conditions. The reader may choose, any combination of vari- 
ables which he considers to be appropriate for the lower strato- 
sphere and develop his own sensitivity table using the results given 
in Fig. 25. 


We showed earlier that E is less than unity because of the 
combined effects of diffusion and mixedness limitations. We now 
examine each of these limitations separately. Figure 31 shows the 
diffusion limitation, Lp . Recall from the definition of 
(Eq. (3^)), that as Lp approaches unity the reaction stops because 
the ambient O 3 cannot be diffused rapidly enough into the plume.' 
Conversely, when Lp approaches 0 , turbulent diffusion supplies 
all the O 3 needed by the reaction. It is clear, in Fig. 3 I, that 
the diffusion limitation is most pronounced for very stable turbu- 
lence fields, very compact plumes and very fast chemical reactions. 


The diffusion' limitation is diminished as 
Increases, or’ k decreases. 


increases, or a. 
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Figure 31. The dependence of the diffusion limitation on 
the turbulence field, source strength, and 
reaction rate constant. The diffusion limita- 
tion Is one component of the efficiency deficit 




Figure 32 shows the effect of k and q 2 on L-j, and clearly 
Illustrates the switch from a diffusion limited reaction for large 
k to a reaction rate-limited condition for small k . 

Figure 33 shows the other deficit to E , the mixedness limita- 
tion, Ljyj . In this figure, the Iso-k curves have a local maximum. 
This peaking phenomenon Is not completely ■ understood at the present 
time. However, the Important point to glean from Fig. 33 Is that 
Ljy[ Is never more than 0.25 For any combination of the Independent 
variables. This suggests that conventional mean value chemistry 
would, at m ost, overpredict the depletion rates for Oo by 25^ by 
neglecting 

Figure 34 shows the sensitivity of Ly, to k for five levels 
of . There appears to be no consistent trends other than that 

L]v[ Is a maximum when the reaction Is neither diffusion rate nor 
reaction rate limited and Ljyj approaches 0 for both limiting 
conditions. A careful comparison of Figs. 32 and 34 shows that 
Lj^ Is always greater than and, therefore, turbulent diffusion 

Is the dominant Influence on the 0 ^ depletion rate. 

We conclude this section with three figures showing E , Lj^ 
and Ljyj contours In the plane. Similar figures were 

previously presented in ref. 2'^but the data obtained during the 
present study permitted an extension of the curves to a larger 
domain and resulted In minor revisions In the original domain. For 
each of these figures, k = 5 xlo 5 (i/pp-sec ) . 

Figure 35 shows E contours as a function of q^ and Cy 
and again demonstrates the low efficiencies associated with low 
turbulence Intensity and compact plumes. For the most part E 
varies more rapidly with Cy than with q2 . 

Figures 36 and 37 show the corresponding Lj^ and Ljyj contours 
respectively. These figures clearly pinpoint the regions of 
maximum diffusion and mixedness limitations which were described 
earlier. 


CONCLUSIONS 


In conclusion, we have demonstrated the sensitivity of the 
coupled chemistry/diffusion model to a wide range of variation of 
the model's Independent variables. We have shown that the effic- 
iency with which the NO^ catalytic cycle destroys ambient Oo Is 
extremely sensitive to the amount of NO emitted and to the relative 
rates of turbulent diffusion and chemical reactions. For repre- 
sentative conditions In the stratosphere, the chemical efficiency 
Is roughly 10 ^. However, a tenfold change In either q 2 ^ k or 
a can produce a factor of five change In the efficiency. If the 
duration of Phase III Is a significant fraction of the total 
residence time of the plume, then these efficiency variations can 
alter 0^ depletion rates by more than a factor of two. 
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Figure 33. The dependence of the tnlxedness limitation on the 
turbulence fields, source strength, and reaction 
rate constant. The mixedness limitation is one 
component of the efficiency deficit. 
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These results point toward those variables which must be 
accurately defined In order to reasonably estimate the Impact of 
SST fleet operations 'on the Oo Inventory. In particular, the need 
for accurate turbulence measurements In the stratosphere Is appar- 
ent. In addition, better estimates of the residence time of the 
plume In the stratosphere which take Into account the aircraft 
flight path (l.e. latitude, longitude and altitude of plume Injec- 
tion) are essential. 

Finally, these results demonstrate the power of the Invariant 
coupled model. By retaining the second-order correlations In the 
controlling equations, we have been able to keep track of the 
separate effects of diffusion, mixedness, and chemistry. We have 
shown that conventional mean-value chemistry can overpredict the 
O 3 depletion rate by at most 25 ^ by neglecting the ralxedness term 

. This inaccuracy may be acceptable and mixedness may be 
neglected but only second-order closure techniques can define the 
magnitude of the error. Even more importantly, we have shown that 
low turbulent diffusion rates are a dominant factor causing low 
efficiency. And^'. therefore, the conventional assumption that the 
reactants are iinlformly mixed at all times Is grossly Inaccurate 
and will significantly overpredict the magnitude of 0 ^ depletion 
due to the SST. ^ 
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